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ABSTRACT 

The  open-literature  material  properties  for  fiber  and  polymeric  matrix,  unit-cell  microstructural  characteristics, 
atomic-level  simulations  and  unit-cell  based  finite-element  analyses  are  all  used  to  construct  a  new  continuum-type 
ballistic  material  model  for  Q]  790  cross-plied  highly-oriented  polyethylene  fiberbased  armor-grade  composite 
laminates.  The  material  model  is  formulated  in  such  a  way  that  it  can 

be  readily  implemented  into  commercial  finite-element  programs  like  ANSYS/Autodyn  [ANSYS/Autodyn 
version  1 1.0,  User  Documentation,  Century  Dynamics  Inc.  a  subsidiary  of  ANSYS  Inc.  (2007)]  and  ABAQUS/Explicit 
[ABAQUS  version  6.7,  User  Documentation,  Dessault  Systems,  2007]  as  a  User  Material  Subroutine.  Model  validation 
included  a  series  of  transient  non-linear  dynamics  simulations  of  the  transverse  impact  of  armor-grade  composite 
laminates  with  two  types  of  projectiles,  which  are  next  compared  with  their  experimental  counterparts.  This  comparison 
revealed  that  a  reasonably  good  agreement  is  obtained  between  the  experimental  and  the  computational  analyses  with 
respect  to:  (a)  the  composite  laminates’  capability,  at  different  areal  densities,  to  defeat  the  bullets  with  different  impact 
velocities;  (b)  post-mortem  spatial  distribution  of  damage  within  the  laminates;  (c)  the  temporal  evolution  of  composite 
armor  laminate  back- face  bulging  and  delamination;  and  (d)  the  existence  of  three  distinct  penetration  stages  (i.e.  an 
initial  filament  shearing/cutting  dominated  stage,  an  intermediate  stage  characterized  by  pronounced  filament/matrix 
de-bonding/decohesion  and  the  final  stage  associated  with  the  extensive  back-face  delamination  and  bulging  of  the 
armor  panel). 
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The  open-literature  material  properties  for  fiber  and  polymeric  matrix,  unit-cell  microstructural  charac¬ 
teristics,  atomic-level  simulations  and  unit-cell  based  finite-element  analyses  are  all  used  to  construct  a 
new  continuum-type  ballistic  material  model  for  0°/90°  cross-plied  highly-oriented  polyethylene  fiber- 
based  armor-grade  composite  laminates.  The  material  model  is  formulated  in  such  a  way  that  it  can 
be  readily  implemented  into  commercial  finite-element  programs  like  ANSYS/Autodyn  [ANSYS/Autodyn 
version  11.0,  User  Documentation,  Century  Dynamics  Inc.  a  subsidiary  of  ANSYS  Inc.  (2007)]  and  ABA- 
QUS/Explicit  [ABAQUS  version  6.7,  User  Documentation,  Dessault  Systems,  2007]  as  a  User  Material  Sub¬ 
routine.  Model  validation  included  a  series  of  transient  non-linear  dynamics  simulations  of  the  transverse 
impact  of  armor-grade  composite  laminates  with  two  types  of  projectiles,  which  are  next  compared  with 
their  experimental  counterparts.  This  comparison  revealed  that  a  reasonably  good  agreement  is  obtained 
between  the  experimental  and  the  computational  analyses  with  respect  to:  (a)  the  composite  laminates’ 
capability,  at  different  areal  densities,  to  defeat  the  bullets  with  different  impact  velocities;  (b)  post-mor¬ 
tem  spatial  distribution  of  damage  within  the  laminates;  (c)  the  temporal  evolution  of  composite  armor 
laminate  back-face  bulging  and  delamination;  and  (d)  the  existence  of  three  distinct  penetration  stages 
(i.e.  an  initial  filament  shearing/cutting  dominated  stage,  an  intermediate  stage  characterized  by  pro¬ 
nounced  filament/matrix  de-bonding/decohesion  and  the  final  stage  associated  with  the  extensive 
back-face  delamination  and  bulging  of  the  armor  panel). 
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1.  Introduction 

Fiber-reinforced  polymer-matrix  composites  are  among  the 
most  advanced  commercially-available  materials.  While  they  are 
widely  used  in  aerospace  and  defense-related  industries,  their 
application  in  construction,  automotive  and  sporting-good  indus¬ 
tries  is  also  quite  common.  The  main  reason  for  the  aforemen¬ 
tioned  widespread  use  of  the  composite  materials  is  their  ability 
to  simultaneously  meet  a  variety  of  functional  and  manufacturing 
requirements.  For  example,  the  new  Boeing  787  Dreamliner  is  pri¬ 
marily  made  of  carbon-fiber-reinforced  epoxy-matrix  composites 
which,  in  addition  to  having  outstanding  mechanical  properties, 
do  not  suffer  from  the  similar  manufacturing  constraints  as  their 
metallic  counterparts/alternatives,  allowing  a  higher  degree  of 
optimization  of  the  787  aerodynamics.  Furthermore,  the  composite 
airframes  weigh  less  and  are  stronger  than  conventional  airframes, 
which  lead  to  improvements  in  the  vehicle’s  operating  efficiency 
and  performance.  Lastly,  carbon-fiber-reinforced  epoxy-matrix 
composites  tend  to  resist  corrosion  and  fatigue,  the  two  phenom¬ 
ena  which  are  well-established  to  cause  gradual  degradation  and 
ultimate  failure  of  metallic  airframes. 
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The  fiber-reinforced  polymer-matrix  composites  like  the  one 
described  above  in  which  the  main  figures  of  merit  are  their  den¬ 
sity-normalized  stiffness  (i.e.  specific  stiffness )  and  density-normal¬ 
ized  strength  (i.e.  specific  strength )  are  commonly  referred  to  as 
“ structural-grade ”  composites.  Many  blast-  and  ballistic-protection 
systems  in  military  and  civilian  applications  are,  on  the  other  hand, 
made  of  another  class  of  fiber-reinforced  polymer-matrix  compos¬ 
ites,  the  so-called  “armor-grade”  composites  [3,4].  The  latter  class 
of  composites  generally  is  optimized  with  respect  to  its  ballistic- 
impact  protection  resistance,  i.e.  with  respect  to  its  energy  absorb¬ 
ing  capability.  Consequently,  the  most  commonly  cited  figures  of 
merit  in  these  materials  are:  (a)  a  critical  level  of  the  projectile’s 
velocity  or  the  projectile’s  kinetic  energy  (generally  referred  to  as 
the  “ ballistic  limit”)  below  which  no  full  perforation  takes  place 
[5,6]  and  (b)  an  extent  to  which  material  ballistic-protection  resis¬ 
tance  is  compromised  in  the  material  systems  which  are  partially 
penetrated  by  projectile(s)  or  whose  strike-face  surface  is  damaged 
by  the  projectile(s). 

The  armor-grade  composites  are  generally  constructed  using 
high  specific-modulus/high  specific-strength  polymeric  fibers  such 
as  aramid  (e.g.  Kevlar®,  Twaron®,  etc.)  or  oriented  polyethylene  fi¬ 
bers  (e.g.  Spectra®,  Dyneema®,  etc.)  with  an  outstanding  impact 
resistance  [7-11].  The  fibers,  in  the  form  of  either  woven  fabrics 
or  in  the  form  of  0°/90°  cross-plied  collimated  continuous 
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filaments,  are  embedded  in  the  resin/polymer  matrix.  To  attain 
maximum  utilization  of  the  inherently-high  transverse-impact 
resistance  of  the  fibers,  the  polymer-matrix  content  does  not 
typically  exceed  20%  by  volume.  As  a  result  of  the  very  low  resin 
content,  these  composites  remain  flexible/compliant  to  relatively 
high  laminate  thicknesses.  Penetration  resistance  of  the  armor- 
grade  composites  is  frequently  increased  through  the  use  of  hybrid 
structures  in  which  a  hard  metallic  or  ceramic  strike-plate  is 
attached  to  the  front  of  an  armor-grade  composite  laminate. 

Armor-grade  composite  laminates  based  on  aramid  fiber-rein- 
forced  phenolic-poly-vinyl-butyral  resin  and  on  0°/90°  cross-plied 
oriented  polyethylene  fiber-reinforced  vinyl-ester  are  widely  used 
in  hard  personnel-armor  systems  (e.g.  protective  helmets)  for  pro¬ 
tection  against  fragments  from  exploding  munitions  [10-15]. 
These  armor-grade  composites  are  also  increasingly  being  used 
for  ballistic  protection  in  light-weight  armored  vehicles,  helicop¬ 
ters,  patrol  boats  and  transportable  shelters  (e.g.  command  shel¬ 
ters)  [10].  Over  the  past  decade,  considerable  efforts  have  been 
invested  in  carrying  out  various  experimental  investigations  in  or¬ 
der  to  identify  and  elucidate  various  penetration-failure  mecha¬ 
nisms  of  the  armor-grade  fiber-reinforced  composites  under 
transverse  impact  loading  and  to  compare  and  contrast  these 
mechanisms  with  those  operating  in  the  related  resin-free  fabrics 
and  resin-rich  structural-grade  composites.  The  main  results  ob¬ 
tained  in  these  investigations  can  be  summarized  as  follows  [lb- 
23]: 

(a)  In  sharp  contrast  to  the  penetration  of  resin-free  fabrics 
which  is  dominated  by  the  successive  fracture  of  individual  yarns 
along  the  periphery  of  the  penetrator  head  and  by  the  side-way/ 
lateral  movement  of  the  yarns  which  enables  them  to  slip  off  from 
the  penetrator,  the  penetration  of  armor-grade  composites  is 
mainly  governed  by  the  failure  of  principal  yarns  (the  yarns  which 
are  in  direct  contact  with  the  penetrator  head).  This  observation  is 
attributed  to  the  effect  of  resin  matrix  on  reducing  yarn  mobility 
which  prevents  them  from  slipping  off  from  the  penetrator.  In  gen¬ 
eral,  stiffer  resin  matrices  (e.g.  vinyl-ester  vs.  polyurethane)  tend  to 
constrain  the  yarn  movement  to  a  greater  degree  and  to  force  the 
penetrator  to  engage  and  fracture  more  yarns  during  penetration. 
This  typically  results  in  improved  ballistic-protection  resistance 
of  armor-grade  composites  and  is  the  reason  that  armor-grade 
composites  reinforced  with  woven-yarn  fabric  are  generally  found 
to  possess  a  higher  energy-absorption  potential  than  their  resin- 
free  fabric  counterparts.  However,  excessive  confinement  of  the 
yarns/filaments  due  to  overly-high  matrix  stiffness  and/or  exces¬ 
sive  amounts  of  the  matrix  may  have  a  deleterious  effect  on  the 
ballistic-protection  performance  of  this  class  of  composites.  The 
latter  effect  is  related  to  the  fact  that  highly  confined  fibers  are 
more  likely  to  fail  in  transverse  shear  before  experiencing  any  sig¬ 
nificant  extensions  in  the  longitudinal  direction. 

(b)  Since  the  energy  absorbed  by  the  armor-grade  composite  is 
found  to  scale  with  the  number  of  broken  yarns  in  its  fabric  con¬ 
stituent,  fiber  tensile  straining  and  ultimate  fracture  is  believed 
to  be  the  dominant  mechanism  for  absorption  of  the  projectile  ki¬ 
netic  energy. 

(c)  In  addition  to  fiber  fracture,  both  woven-fabric-reinforced 
and  cross-plied  fiber-reinforced  composite  laminates  are  generally 
found  to  include  additional  complex  failure  processes  such  as:  (i) 
delamination,  (ii)  a  plug  punch  out,  (iii)  resin  matrix  cracking, 
and  (iv)  fiber  pull-out.  These  failure  modes  are  also  typically  ob¬ 
served  in  conventional  structural-grade  composites  reinforced 
with  glass  or  carbon  fibers. 

(d)  In  the  case  of  multi-ply  armor-grade  composite  laminates 
reinforced  with  either  cross-plied  collimated  Spectra  fibers  or  with 
woven  Spectra  fabrics,  the  following  fracture  modes  are  most  often 
observed  [1]:  (i)  sequential  delamination,  (ii)  plug  punch  out  in¬ 
duced  by  the  through-the-thickness  shear,  and  (iii)  combined  fiber 


shearing/cutting  and  fiber  tensile  failure.  In  the  cross-plied  lami¬ 
nates,  fibers  in  the  top  plies  are  typically  found  to  fail  by  shear¬ 
ing/cutting,  primarily  along  the  edges  of  the  projectile.  Fibers 
located  in  the  back  layers  of  the  laminates,  on  the  other  hand,  gen¬ 
erally  fail  in  tension  although  in  thin  laminates,  the  lateral  motion 
of  fibers  and/or  fiber  pull-out  rather  than  fiber  tensile  straining  to 
fracture  is  sometimes  observed. 

(e)  The  delamination  in  the  cross-plied  Spectra  fiber-reinforced 
composite  laminates  is  typically  found  to  resemble  the  “generator 
strip ”  phenomenon  [22]  seen  in  glass  fiber-reinforced  epoxy-ma¬ 
trix  structural-grade  composites.  That  is,  under  transverse  impact, 
the  projectile  pushes  a  strip  of  the  first  lamina  toward  the  rear  of 
the  laminate  which  induces  shear  cracks  in  the  resin  matrix  paral¬ 
lel  to  the  fibers  and  applies  a  transverse  load  to  the  second  lamina. 
This,  in  turn,  causes  separation  between  the  first  two  laminae,  i.e. 
delamination.  After  the  aforementioned  delamination  process  had 
taken  place  successively  through  the  entire  thickness  of  the  lami¬ 
nate  via  the  same  mechanism  and  penetration  of  the  laminate 
has  occurred,  narrow  strips  of  damage  zone  remain  visible  under 
transmitted  light  and  the  strips  are  found  to  tend  to  follow  the 
respective  fiber  orientation  in  the  panel.  These  strips  typically  con¬ 
tain  numerous  matrix/fiber  interface  cracks.  In  addition,  a  circular 
delamination  zone  is  generally  seen  around  the  perforation  hole. 

(f)  In  contrast  to  the  case  of  cross-plied  fiber-reinforced  com¬ 
posite  laminates,  fabric-reinforced  laminates  are  found  to  exhibit 
much  less  lateral  movements  of  reinforcing  fibers  during  the  pen¬ 
etration  of  the  projectile  [3,4].  Even  in  thin  panels,  fibers  appar¬ 
ently  failed  due  to  shearing/cutting  in  the  laminae  near  the 
strike-face  and  in  tension  at  the  rear  of  the  completely  penetrated 
laminates.  The  presence  of  a  narrow  strip  of  the  first  lamina  pushed 
forward  by  the  penetrator  is  generally  not  observed.  Instead,  the 
delamination  zones  are  observed  preferentially  along  the  two  rein¬ 
forcement  directions  of  woven  fabric.  However,  these  damage 
zones  are  closely  integrated  with  the  circular  delamination  zone 
around  the  perforation  hole.  The  occurrence  of  less  anisotropic  pat¬ 
tern  of  delamination  was  linked  with  the  presence  of  resin-rich 
pockets  between  the  reinforcing  layers  and  with  a  greater  con¬ 
straint  to  matrix  crack  propagation  parallel  to  the  fibers/yarns. 

(g)  Up  to  the  thickness  of  ~  3  mm,  the  kinetic  energy  for  full 
perforation  of  armor-grade  composites  has  been  found  to  depend 
on  the  laminate  thickness  in  a  way  similar  to  that  observed  in  duc¬ 
tile  monolithic  materials  (e.g.  poly-carbonate  or  aluminum). 

The  full-perforation  kinetic  energy  vs.  laminate  thickness  rela¬ 
tionship,  however,  was  found  to  be  somewhat  non-linear.  This 
finding  has  been  attributed  to  the  unique  mode  of  tensile  failure 
seen  in  these  materials  for  which  the  critical  level  of  kinetic  energy 
for  full  perforation  is  lowered  by  the  fiber/yarn  mobility. 

The  first  use  of  fiber-based  composites  (primarily  nylon  ( poly¬ 
amide )  fabric  and  E-glass  fiber/ethyl  cellulose  composites)  in  body 
armor  systems  in  place  of  the  traditionally  used  metallic  solutions 
can  be  traced  back  to  the  Korean  War  [24].  Although,  primarily  due 
to  their  low  cost,  nylon  and  E-glass  fibers  are  still  being  used  today, 
high-performance  polymeric  fibers  are  now  the  standard  in  most 
fiber-reinforced  body-armor  applications.  The  high-performance 
polymeric  fibers  used  today  are  characterized  by  substantially  im¬ 
proved  strength,  stiffness  and  energy-absorbing  capacity.  Among 
these  high-performance  fibers  the  most  notable  are:  (a)  poly-ara- 
mids  (e.g.  Kevlar®,  Twaron®,  Technora®);  (b)  highly-oriented  ultra 
high  molecular  weight  polyethylene ,  UHMWPE  (e.g.  Spectra®, 
Dyneema®);  (c)  poly-benzobis-oxazole ,  PBO  (e.g.  Zylon®);  and  (d) 
poly-pyridobisimi-dazole ,  PIPD  (e.g.  M5®).  When  tested  in  tension, 
all  these  materials  differ  significantly  from  the  nylon  fibers,  having 
very  high  absolute  stiffness,  extremely  high  specific  strength,  and 
quite  low  (<4%)  strains-to-failure.  These  fibers  essentially  behave, 
in  tension,  as  rate-independent  linear  elastic  materials.  When 
tested  in  transverse  compression,  however,  these  fibers  are  similar 
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to  nylon  and  can  undergo  large  plastic  deformation  without  a  sig¬ 
nificant  loss  in  their  tensile  load-carrying  capacity.  This  behavior  is 
quite  different  from  that  found  in  carbon  or  glass  fibers,  which  tend 
to  shatter  under  transverse  compression  loading  conditions. 

The  ballistic  performance  of  high-performance  polymeric  fibers 
is,  in  general,  quantified  with  respect  to  their  ability  to:  (a)  absorb 
the  projectile’s  kinetic  energy  locally  and  (b)  spread  out  the  ab¬ 
sorbed  energy  fast  before  local  conditions  for  the  failure  are  met. 
In  simple  terms,  the  ability  of  high-performance  fibers  to  absorb 
energy  per  their  unit  mass,  ESP,  is  related  to  the  fiber  tenacity/fail¬ 
ure  strength,  the  fiber  strain-to-failure,  £/ai/,  and  the  fiber  den¬ 
sity,  p ,  as: 

Esp  —  0  ■  5  O' fan  Bfail  /  P  (1) 

The  ability  of  fibers  to  spread  out  energy  is  governed  by  their 
speed  of  sound,  vS0Und ,  which  is  defined  in  terms  of  their  axial  mod¬ 
ulus  of  elasticity,  E,  and  their  density  as: 

Vsound  =  (E/pf/2)  (2) 

In  Fig.  1,  the  two  aforementioned  ballistic  performance  param¬ 
eters  are  displayed  for  the  most-commonly  used  high-performance 
fibers.  A  summary  of  the  key  properties  of  the  same  set  of  high- 
performance  fibers  is  provided  in  Table  1. 

While  the  results  displayed  in  Fig.  1  clearly  reveal  a  high  bal¬ 
listic  potential  of  the  high-performance  fibers  in  general  (and 
specifically  of  the  highly-oriented  UHMWPE  fibers,  the  type  of  fi¬ 
ber-reinforcements  considered  in  the  present  work),  full  utiliza¬ 
tion  of  this  potential  in  armor-grade  composites  turned  out  to 
be  a  formidable  challenge  because  a  number  of  additional  factors 
(e.g.  fabric/ply  structure/architecture,  ply  areal  density,  fiber-to- 
fiber/yarn-to-yarn  and  fiber/yarn-to-projectile  friction,  type  of 
polymeric  matrix,  composite  processing  and  fabrication  condi¬ 
tions,  shape,  mass  and  mechanical  properties  of  the  projectile 
to  be  defeated,  etc.)  become  important.  To  overcome  these  chal¬ 
lenges,  the  development  of  flexible-armor  systems  has  started  to 
rely  increasingly  more  on  the  use  of  transient  non-linear  dynam¬ 
ics  computational  analyses  of  the  ballistic  response  of  armor 
when  impacted  with  high-speed  projectiles.  For  these  analyses 
to  yield  reliable  predictions  and  for  them  to  be  used  as  comple¬ 
ments  to  the  accompanying  experimental  investigations,  high- 
fidelity  physically-based  material  models  for  the  armor-grade 
composite  materials  must  be  available. 


Mass-basedEnergyAbsorptionCapacity,kJ/kg 


Fig.  1.  Sound  speed  vs.  mass-based  energy  absorption  capacity  for  a  number  of 
high-performance  fibers. 


A  review  of  the  public-domain  literature  carried  out  as  part  of 
the  present  work  revealed  the  existence  of  several  material  models 
for  armor-grade  composites  [25-30].  While  such  models  have  pro¬ 
vided  an  important  insight  into  the  roles  of  a  number  of  factors 
mentioned  above,  they  suffer  from  three  major  shortcomings:  (a) 
some  are  more  phenomenological,  i.e.  less  physically-based  in 
their  character;  (b)  others  require  the  knowledge  of  a  relatively 
large  number  of  parameters;  and  (c)  may  not  very  efficient  compu¬ 
tationally.  In  addition  to  the  models  mentioned  above,  purely  phe¬ 
nomenological  models  [e.g.  31]  also  exist  in  the  literature.  Such 
models  are  the  result  of  extensive  experimental  efforts  and  typi¬ 
cally  have,  within  the  same  family  of  armor-grade  composite 
materials,  a  high  practical  utility.  However,  they  provide  no  insight 
into  the  complicated  physics  of  projectile/armor  interactions  and 
can  not  be  used  across  the  boundaries  of  different  armor-type  com¬ 
posite  families. 

To  overcome  the  aforementioned  limitations  of  the  two  groups 
of  material  models,  a  new  multi-scale  physically-based  computa- 
tionally-efficient  material  model  for  UHMWPE-filament  (e.g.  Spec¬ 
tra®,  Dyneema®,  etc.)  based  armor-grade  composites  is  developed, 
parameterized  and  validated  in  the  present  work.  A  preliminary 
version  of  this  model  without  details  regarding  the  atomic-level 
analysis  of  the  filament/matrix  interfacial  bonding/de-bonding 
and  without  detailed  validation  was  previously  reported  Grujicic 
et  al.  [32].  Since  it  is  generally-established  that  for  the  UHMWPE 
fiber-based  armor-grade  composites,  a  substantially  higher  ballis¬ 
tic-protection  performance  is  obtained  when  such  fibers  are  used 
as  0°/90°  cross-plied  layers  of  highly-oriented  filaments  rather 
than  woven  fabrics,  only  the  former  composite-laminate  architec¬ 
ture  will  be  addressed  in  the  present  work.  In  passing,  it  should  be 
mentioned  that  it  is  believed  that  the  deflection  of  stress  waves  at 
the  yarn/yarn  or  fiber/fiber  cross-ever  points  in  woven  fabric  (the 
process  which  lowers  the  ability  of  fibers  to  spread  out  energy 
along  their  axis)  is  the  main  reason  for  their  inferior  ballistic 
performance. 

As  stated  above,  the  main  objective  of  the  present  work  was  to 
develop,  parameterize  and  validate,  (against  the  relevant  experi¬ 
mental  results),  a  simple  physically-based  computationally-effi- 
cient  continuum-level  material  model  for  a  prototypical  0°/90° 
cross-plied  oriented  polyethylene  fiber-based  armor-grade  com¬ 
posite  material.  While  developing  this  model,  unit-cell  level  fi¬ 
nite-element  analyses  of  the  meso-scale  material  mechanical 
response  and  properties  and  atomic-level  simulations  of  the  fila¬ 
ment/matrix  bonding/de-bonding  had  to  be  employed  giving  the 
present  approach  a  multi-length  scale  character. 

It  should  be  noted  that  within  a  fully  multi-scale  computational 
approach,  the  underlying  boundary-value  problem  is  set-up  and 
simultaneously  solved  at  several  length  scales  (e.g.  at  atomistic, 


Table  1 

Typical  mechanical  properties  of  high-performance  fibers. 


Fiber 

type 

Failure  strength 
(GPa) 

Failure 

strain 

Axial  modulus 
(GPa) 

Density 

(kg/m3) 

Aramid 

2. 8-3.2 

0.015- 

0.045 

60-115 

1390-1440 

HMWPE 

2. 8-4.0 

0.029- 

0.038 

90-140 

970-980 

LCP 

2. 7-2.9 

0.033- 

0.035 

64-66 

1400-1420 

PBO 

5.4-5. 6 

0.024- 

0.026 

270-290 

1540-1560 

PI  PD 

3.9-4.1 

0.011- 

0.013 

320-340 

1690-1710 

Nylon 

0.06-0.08 

1. 5-2.5 

1.0-1. 5 

1070-1170 

S-glass 

4.64-4.66 

0.053- 

0.055 

82-92 

2470-2490 
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meso,  macro).  This  approach  will  be  adopted  in  our  future  work.  In 
the  present  investigations,  a  simpler,  and  computationally  more 
efficient  yet  less  accurate  approach  was  adopted.  That  is,  results 
from  the  experimental  investigation  of  Iremonger  [36]  are  used 
to  infer  the  main  deformation  and  fracture  mechanisms  of  the 
composite  material  under  investigation  and  a  coupled  meso- 
scale/continuum  level  model  is  constructed.  However,  critical  fila¬ 
ment/matrix  de-bonding  parameters  needed  in  this  model  were 
determined  in  a  separate  atomic  level  computational  investigation, 
described  in  Appendix  A. 

The  organization  of  the  paper  is  as  follows:  Details  regarding 
the  computational  procedures  employed  to  develop  a  new  unit-cell 
continuum-damage  based  material  model  for  a  prototypical  0°/90° 
cross-plied  uni-directional  UHMWPE-filament  based  armor-grade 
composite  and  the  implementation  of  this  model  into  a  material 
user  subroutine  suitable  for  use  in  commercial  finite-element  pro¬ 
grams  are  presented  in  Section  2.  The  formulation  of  a  simple  pro¬ 
jectile-target  impact  problem  used  to  validate  the  new  material 
model  is  described  in  Section  3.  Main  results  obtained  in  the  cur¬ 
rent  work  are  presented  and  discussed  in  Section  4.  The  main  sum¬ 
mary  points  and  conclusions  resulting  from  the  present  work  are 
listed  in  Section  5.  A  brief  discussion  regarding  the  atomic-scale 
computational  analysis  used  to  investigate  filament/matrix  bond¬ 
ing/de-bonding  behavior  is  presented  in  Appendix  A. 

2.  Material  model  formulation  and  implementation 

In  this  section  and  its  subsections,  a  detailed  account  is  given  of 
the  procedure  used  to  develop  a  new  unit-cell  continuum-damage 


Fig.  2.  The  relationship  between  a  unit  cell  and  the  corresponding  material  point  in 
an  anisotropic  continuum. 


based  multi-scale  material  model  for  a  prototypical  single-lamina 
0°/90°  cross-plied  uni-directional  UHMWPE-filament  based  ar¬ 
mor-grade  composite.  Also  details  regarding  the  implementation 
of  the  model  into  a  material-user  subroutine  suitable  for  use  in 
commercial  finite-element  packages  are  presented.  The  basic  idea 
behind  the  unit-cell  based  approach  is  that  the  mechanical  re¬ 
sponse  of  the  unit-cell  (consisting  of  high-stiffness/high-strength 
polymeric  filament  segments  and  a  compliant  polymeric  matrix) 
can  be  smeared  out  (homogenized)  into  an  equivalent  response 
of  a  (anisotropic)  continuum  material.  A  simple  schematic  of  the 
unit  cell  which  is  used  to  represent  0°/90°  cross-plied  unidirec¬ 
tional  UHMWPE-filament  based  armor-grade  composites  allotted 
to  a  single  filament  crossover  is  depicted  in  Fig.  2(a).  Its  contin¬ 
uum-level  material  point  counterpart  is  represented  in  Fig.  2(b). 
Within  the  continuum-material  framework,  filaments  are  not  rep¬ 
resented  explicitly  but  rather  by  two  material  directions  whose 
orientations  are  denoted  in  terms  of  material  vectors,  g1  and  g2. 
(Please  note  that  vectors  are  denoted  using  a  bold  lower-case  font, 
tensors  using  a  bold  upper-case  font  while  scalars  using  a  non-bold 
font.)  The  “unit-cell”  term  is  used  to  denote  the  basic  structural 
block  so  that  a  piece  of  the  armor-grade  composite  material  can 


Fig.  3.  Typical  finite-element  meshes  used  in  the  unit-cell  computational  analyses 
to  discretize:  (a)  the  two  filaments  and  (b)  the  matrix. 


472 


M.  Grujicic  et  al./  Composites:  Part  B  40  (2009)  468-482 


Table  2 

The  orthotropic  linear  elastic  material  data  for  UHMWPE  filaments  [31]. 


En 

E22 

E33 

G\2 

G13 

G23 

Vl2  V13  V23 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

118.0 

6.0 

6.0 

6.0 

6.0 

6.0 

0.3  0.3  0.4 

Axial  failure  strain  =  0.05;  transverse  shear  strength  =  350  MPa. 


be  considered  as  a  result  of  the  repetition  of  this  block  in  three 
orthogonal  directions. 

Coupling  between  the  continuum  material  formulation  and  the 
unit-cell  geometry  and  mechanical  response  is  done  in  the  follow¬ 
ing  way:  (a)  the  deformation  state  of  a  continuum  material  point 
(as  quantified  by  the  corresponding  deformation  gradient)  is  used 
to  update  the  unit-cell  geometry;  (b)  the  updated  unit-cell  geom¬ 
etry  and  the  state  of  the  continuum  material  at  the  end  of  the  pre¬ 
vious  time  increment  are  used  to  update  the  extent  of  structural 
damage  in  the  unit  cell;  and  (c)  the  updated  material  state  ob¬ 
tained  in  point  (b)  is  then  used  to  compute  the  stress  state  at  the 
end  of  the  current  time  increment. 

It  must  be  recognized  that  in  order  for  the  aforementioned 
approach  to  be  valid  (i.e.  in  order  for  homogenization  of  the  ar¬ 
mor-grade  composite  unit-cell  response  to  be  justified),  the 
characteristic  lengths  in  the  numerical  analysis  in  which  the  model 
is  used  (e.g.  the  projectile  and  the  target  dimensions  in  a  projectile/ 
target  impact  problem  analyzed  in  the  present  work  and  the  asso¬ 
ciated  stress/strain  gradient  ranges)  must  be  large  in  comparison 
to  the  unit-cell  edge  lengths.  In  most  practical  situations,  this  ap¬ 
pears  to  be  the  case  since  the  unit  cell  edge  length  is  in  the  range 
between  10  and  30  pm. 

2.1.  Meso-scale  unit-cell  level  finite -element  analyses 

The  salient  feature  of  the  proposed  computational  approach  is 
that  the  mechanical  response  of  a  continuum-level  material  point 
(corresponding  to  a  unit  cell  in  the  armor-grade  composite)  and 
the  accompanying  changes  in  constituent  materials  (primarily 
those  associated  with  the  filament/matrix  interfacial  de-bonding) 
can  be  inferred  by  carrying  out  a  series  of  meso-scale  finite-ele¬ 
ment  analyses  pertaining  to  relatively  simple  mechanical  tests 
of  the  unit  cell.  In  these  analyses,  a  detailed  representation  of 
the  unit-cell  microstructure  is  considered.  In  this  section,  details 
are  presented  regarding  the  geometrical  models  used  in  the  con¬ 
struction  of  the  unit  cell,  material  properties/models  assigned  to 
the  filament  segments,  matrix  and  the  filament/matrix  interfacial 
bonding/de-bonding,  and  the  finite-element  analyses  used  to 
determine  the  mechanical  response  and  the  material  evolution 
under  different  loading  conditions.  As  will  also  be  shown  in  this 
section,  to  obtain  quantitative  information  about  the  filament/ma¬ 
trix  interfacial  bonding/de-bonding,  atomic-level  simulations 
were  employed.  This  portion  of  the  work  is  covered  in  more  de¬ 
tails  in  Appendix  A.  It  should  be  noted  that,  as  correctly 
pointed  out  by  one  of  the  reviewers  of  the  present  work,  interfa¬ 
cial  de-bonding  parameters  can  be  also  measured  experimentally 
using  one  of  the  methods  such  as  the  so-called  “fiber  push-out” 
method. 


An  example  of  the  finite-element  meshes  used  in  the  unit-cell 
computational  analyses  is  displayed  in  Fig.  3(a)  and  (b).  20,847  first 
order  tetrahedron  elements  (ABAQUS/Explicit  designation  C3D4) 
are  used  to  discretize  each  of  the  two  filament  segments, 
Fig.  3(a),  while  21,606  elements  of  the  same  type  are  used  to  dis¬ 
cretize  the  matrix,  Fig.  3(b).  Bonding  between  the  matrix  and  the 
filaments  is  represented  using  7056  “ cohesive ”  elements,  (ABA¬ 
QUS/Explicit  designation  COH3D6). 

The  polymeric  filaments  (assumed  to  be  based  on  the 
UHMWPE)  are  modeled  as  orthotropic  (more  precisely  as  planar 
isotropic)  linear  elastic  materials  (up  to  the  point  of  failure  under 
axial  tension  or  transverse  shear)  with  the  unique  material  direc¬ 
tion  being  aligned  with  the  filament  axis.  A  summary  of  the  elastic 
and  failure  properties  of  the  filament  material  is  provided  in  Table 
2. 

The  polymeric  matrix  (assumed  to  be  based  on  styrene-iso- 
prene-styrene  tri-block  copolymer  [33])  is  modeled,  due  to  atten¬ 
dant  high-deformation  rate  conditions,  as  a  linear  isotropic 
material  with  a  Young’s  modulus  of  3  GPa  and  a  Poisson’s  ratio  of 
0.4. 

Bonding  between  the  filaments  and  the  matrix  is  modeled  using 
traction  vs.  interfacial  displacement-discontinuity  relations  (one 
for  the  normal  and  one  for  the  tangential  displacements).  These 
two  relations  are  characterized  by  a  linear  traction  vs.  displace¬ 
ment/discontinuity  relation  unto  the  point  of  damage  initiation 
and  with  a  linear  “ Downhill ”  post  damage  relationship.  Conse¬ 
quently  the  two  modes  “Normal  and  Shear ”  of  interfacial-bonding 
damage  are  each  characterized  by  three  parameters:  (a)  critical 
normal  or  shear  interfacial-displacement  discontinuities  at  which 
damage  initiation  begins;  (b)  the  corresponding  normal  or  shear 
interfacial  strengths;  and  (c)  normal  or  shear  interfacial-displace- 
ment  discontinuities  at  which  complete  filament/matrix  decohe¬ 
sion  takes  place.  A  summary  of  the  interfacial  cohesion 
parameters  used  in  the  present  work  is  given  in  Table  3.  These 
parameters  were  determined  in  a  separate  molecular-statics  based 
investigation  of  atomic-level  mechanical  properties  of  the  compos¬ 
ite  materials  consisting  of  the  unidirectional  UHMWPE  filaments 
and  an  amorphous  polymeric  matrix.  Few  details  regarding  the 
atomic-level  analysis  carried  out  and  the  results  obtained  are  pre¬ 
sented  in  Appendix  A.  The  procedure  used  closely  follows  that  pre¬ 
sented  in  our  recent  work  [34]  in  which  atomic  level  properties  of 
composite  materials  consisting  of  multi-walled  carbon  nanotube 
reinforcements  and  a  poly-vinyl-ester-epoxy  matrix  were 
investigated. 

Interactions  between  the  filaments  and  the  matrix  after  decohe¬ 
sion  are  accounted  for  through  a  “Hard"  pressure  vs.  over  closure 
algorithm  within  which  the  interacting  bodies  must  be  in  contact 
before  they  can  interact  and  the  pressure  levels  that  can  be  trans¬ 
mitted  through  the  contact  interactions  are  unbounded.  Relative 
sliding  of  the  contacting  bodies  is  opposed  by  a  frictional  force 
based  on  a  constant  friction  coefficient. 

The  following  simple  mechanical  tests  were  carried  out  using 
the  meso-scale  unit-cell  based  finite-element  approach  described 
above:  (a)  uniaxial  tension  along  the  axis  of  one  of  the  filaments 
(i.e.  along  the  directions  1  or  2);  (b)  uniaxial  tension  in  a  direction 
normal  to  the  single-lamina  surface  (direction  3);  (c)  in-plane  1-2 
shear;  and  (d)  the  transverse  shear. 


Table  3 

Normal  and  shear  filament/matrix  de-bonding  parameters  used  in  the  present  work. 


Normal  de-bonding 

Shear  de-bonding 

Initiation  displacement- 
discontinuity  (pm) 

Bond 

strength 

(MPa) 

Complete  de-bonding  displacement- 
discontinuity  (pm) 

Initiation  displacement- 
discontinuity  (pm) 

Bond 

strength 

(MPa) 

Complete  de-bonding  displacement- 
discontinuity  (pm) 

0.05 

18.0 

1.1 

0.9 

23.0 

2.1 
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It  is  well-established  that  the  presence  of  free  surfaces  (along 
which  loading  is  applied)  in  isolated  unit  cells,  like  the  one  used 
in  the  present  work,  can  alter  the  unit  cell  mechanical  response  rel¬ 
ative  to  the  one  displayed  by  “bulk”  unit  cells  subjected  to  the 
same  mode  of  loading.  To  provide  the  first-order  assessment  of 
these  effects,  few  mechanical-test  simulations  were  carried  out 
using  a  block  of  3  x  3  x  3  =  27  unit  cells.  The  results  obtained  show 
that  the  “surface-loading  effects”  can  be  significant.  However, 
these  effects  were  not  included  in  the  current  rendition  of  the 
material  model  since  this  would  have  required  the  level  of  effort 
not  available  in  the  present  work. 

For  each  of  the  aforementioned  tests,  a  series  of  loading- 
unloading-reloading  cycles  was  applied  in  order  to  detect  the  on¬ 
set  of  interfacial  de-bonding  and  the  resulting  degradation  in  the 
corresponding  continuum-material  stiffness  parameters.  To  deter¬ 
mine  how  one  mode  of  loading  may  affect  all  the  unit-cell  stiffness 
parameters,  loading  is  done  in  one  mode  while  subsequent  reload¬ 
ing  is  done  in  all  the  modes  (one  at  a  time).  It  should  be  noted  that 
simulations  of  the  mechanical  test  mentioned  above  were  done 
under  loading  rates  comparable  to  those  typically  encountered  in 
the  bullet/target  impact  problems.  Consequently,  the  inertial  ef¬ 
fects  were  accounted  for,  at  least  to  the  first  order  of 
approximation. 

The  results  obtained  suggest  that  interfacial  de-bonding  is 
mainly  caused  by  the  through-the-thickness  tension  (in  direction 
3)  and  by  the  in-plane  shear  (1-2  shear)  and  that  F33,  Gu,  G23 
and  G3i  are  mostly  degraded  by  interfacial  de-bonding.  In  addition, 
these  four  stiffness  parameters  are  found  to  degrade  essentially 
linearly  with  the  extent  of  interfacial-bonding  damage,  D. 

In  accordance  with  the  aforementioned  observations  and  the 
atomic-level  interfacial  de-bonding  initiation  parameters  listed  in 
Table  3,  the  following  strain-based  damage  initiation  criterion 
was  derived: 


Fig.  4.  A  schematic  of  the  computational  procedure  used  to  determine  the  onset 
and  track  the  progress  of  material  damage  (interfacial  de-bonding)  within  the  unit 
cell. 


present  model  is  not  suitable  for  use  in  the  calculations  dealing 
with  low  loading-rate  structural  response  of  these  materials. 

2.2.  Determination  of  the  unit-cell  current  geometry  and  architecture 


-^-Y+ 

f33,initj 


1 
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where  e33,imt  and  yuinit  are  pure  normal  and  shear  strains  at  which 
damage  initiation  is  first  observed.  In  a  vs.  -112-  plot,  Eq.  (3)  de- 
fines  a  unit  {failure -initiation)  circle.  Within  the  same  plane,  the  con¬ 
dition  at  which  complete  damage-induced  degradation  takes  place 
is  defined  by  an  ellipse  in  the  form: 


(4) 


where  s33jaii  and  y12/ai7  are  pure  normal  and  shear  strains  at  which 
complete  degradation  (D  =  1 )  takes  place. 

When  the  (£33,712)  strain  state  of  a  material  point  (e.g.  point  B 
in  Fig.  4)  lies  between  the  damage  initiation  circle  and  the  failure 
ellipse,  the  corresponding  extent  of  material  damage  is  defined 
as  the  ratio  of  line  segments  AB  and  AC  indicated  in  Fig.  4.  The  four 
continuum  level  damage  parameters  s33jinit,  s33fau,  yu,mt  and  7i2m 
are  determined  using  the  unit-cell  finite-element  analyses  de¬ 
scribed  in  this  section.  The  values  of  these  four  parameters  are: 
0.01,  0.04,  0.01  and  0.04,  respectively. 

It  should  be  noted  that,  as  correctly  pointed  out  by  one  of  the 
reviewers  of  the  present  work,  the  mechanical  response  (including 
the  failure  behavior)  of  composite  materials  is  loading-rate  depen¬ 
dent.  Since  the  present  material  model  is  intended  to  be  used  only 
in  the  high  loading-rate  regime  (the  regime  which  dominates  bul¬ 
let/target  interaction  conditions),  under  which  the  constituent 
materials  tend  to  behave  essentially  as  linear  elastic  with  weak 
rate  dependency,  the  material  model  parameters  utilized  in  the 
present  work  are  treated  as  rate  independent.  Consequently,  the 


As  discussed  earlier,  a  critical  step  in  the  development  of  the 
present  continuum-damage  material  model  is  establishment  of 
the  relationship  between  the  continuum-level  material-point 
deformation  state  and  the  meso-scale  unit-cell  geometry.  This  to¬ 
pic  is  covered  in  the  present  section. 

In  general  six  independent  geometrical  parameters  are  needed 
to  fully  describe  the  current  geometry  of  the  unit  cell.  These 
parameters  include:  (a)  The  three  unit-cell  edge  lengths, 
a,  (i  =  1—3);  (b)  the  in-plane  shear  inter-filament  included  angle, 
6 ;  and  (c)  the  two  out-of-plane  shear  angles,  and  i/r.  In  this  section 
it  is  shown  how  these  parameters  are  related  to  the  continuum-le¬ 
vel  deformation  state  of  the  material  point  corresponding  to  the 
unit  cell  in  question. 

At  the  continuum  level,  the  state  of  deformation  at  a  given 
material  point  is  described  by  the  deformation  gradient,  F,  whose 
components  in  a  Cartesian  coordinate  system  are  defined  in  Eq. 
(5)  in  Ref.  [32]:  Next,  also  at  the  continuum  level,  the  0°  and  90° 
filaments  can  be  described  using  vectors  a,  (i  =  1 , 2)  aligned  with 


Table  4 

Experimental  [36]  and  the  corresponding  computational  [Present  Work]  results 
pertaining  to  the  success  of  armor-grade  composite  test  panel  to  stop  a  M855  bullet  at 
different  initial  bullet  velocities. 


Test-panel 

Areal  density 

Bullet  velocity  (m/s) 

thickness  (mm) 

(kg/m2) 

600 

700 

800 

900 

4.2 

4 

- 

- 

- 

G/G 

11 

10.5 

G/G 

G/G 

G/G 

G/G 

15 

13.7 

- 

O/O 

O/O 

G/G 

22 

21 

U/U 

U/U 

U/U 

U/U 

32 

31 

u/u 

U/U 

U/U 

U/U 

Nomenclature:  U  -  undermatched,  O  -  overMatched,  G  -  grossly  overmatched; 
Experiment/Computation. 
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the  axis  of  these  filaments  and  the  length  of  these  vectors  can  be 
set  equal  to  the  corresponding  current  in-plane  unit-cell  edge 
lengths,  a,  (i=  1,2).  These  vectors  and  the  vector  a3  which  is 
aligned  with  the  out-of-plane  unit-cell  edge  can  be  related  to  their 
initial  counterparts,  ai0  (i  =  1—3)  as  defined  in  Eq.  (6)  in  Ref.  [32]. 
The  length  of  each  a*  (i  =  1—3)  can  be  defined  using  Eq.  (7)  in  Ref. 
[32]  while  the  inter-filament  included  angle  6  can  be  computed 
using  Eq.  (7)  in  Ref.  [32]:  An  equation  analogous  to  Eq.  (8)  in  Ref. 
[32]  can  next  be  used  to  define  the  out-of-plane  shear  angles  </> 
and  i j/.  Eqs.  (5)-(8)  in  Ref.  [32]  can  be  used  to  show  that  the  current 
geometry  of  the  unit  cell  can  be  determined  using  the  original  unit¬ 
cell  edge  lengths  and  the  current  value  of  the  deformation 
gradient. 

Once  the  current  unit-cell  parameters  are  defined,  standard 
relations  are  invoked  to  compute  the  corresponding  normal  and 
shear  strains.  Next,  through-the-thickness  normal  strain  s33  and 
the  in-plane  shear  strain  y12  and  the  procedure  outlined  in  the  pre¬ 
vious  section  are  used  to  update  the  extent  of  material  damage  and 
the  affected  stiffness  moduli.  It  should  be  noted  that  material  dam¬ 
age  is  irreversible,  i.e.  D  cannot  decrease  during  the  deformation 
history  of  a  material  point. 

2.3.  Determination  of  the  continuum  material-point  stress  state 

Once  the  extent  of  material  damage  is  updated  then 
E33,  G12,  G23  and  G3i  are  degraded  by  multiplying  their  initial  val¬ 
ues  by  a  factor  (1.0-cD),  where  0.0  ^  c  ^  1.0  is  an  elastic-modulus 
dependent  parameter.  Since  the  continuum  material  is  modeled  as 
a  linear  elastic  orthotropic  material  with  degradable  stiffness  mod¬ 
uli,  the  standard  relationships  are  used  to  compute  the  stress  com¬ 
ponents  from  the  updated  strain  components  and  the  updated 
material  stiffness  matrix. 

Once  the  stresses  are  updated,  the  occurrence  of  filament  failure 
is  investigated.  Filaments  are  allowed  to  fail  in  one  of  the  two  fol¬ 
lowing  modes:  (a)  in  tension,  when  the  tensile-strain/tensile-stress 
reaches  a  critical  value  or  (b)  due  to  transverse  shear,  when  the 
corresponding  transverse  shear  stress  reaches  a  critical  value. 
When  either  of  these  two  filament  failure  modes  takes  place,  the 
corresponding  in-plane  normal  stress(es)  and  the  corresponding 
transverse  shear  stress  are  set  to  a  small  residual  value  associated 
with  the  remaining  matrix  ligaments.  Once  the  stresses  are  up¬ 
dated  to  include  the  effect  of  filament  failure  they  are  ready  to 
be  returned  to  the  finite-element  solver  for  the  computation  of 
the  global  equilibrium. 

The  present  material  model  is,  thus,  constructed  in  such  a  way 
so  that  it  can  account  for  the  competition  between  the  following 
two  processes:  (a)  Transverse-shear  loading  which  is  promoted 


by  good  filament/matrix  bonding  and  higher  matrix  stiffness.  If 
sufficiently  high  transverse  shear  stresses  are  developed  they  can 
cause  shear/type  failure  of  the  filament(s).  In  this  case,  the  energy 
absorbed  by  the  filaments  is  relatively  small  and,  consequently, 
ballistic/protection  performance  of  the  armor-grade  composite 
laminate  is  inferior  and  (b)  stretching  of  the  filaments  till  the  point 
of  failure.  This  process  is  promoted  by  filament/matrix  de-bonding 
which  enables  the  filaments  to  deform  independently  of  the  ma¬ 
trix.  In  this  case,  the  energy  absorbed  by  the  filaments  is  maximum 
and  the  ballistic/protection  of  the  armor  is  greatly  enhanced.  It 
should  be  noted  that  some  critical  level  of  filament-matrix  bonding 
is  needed  to  ensure  that  the  filaments  are  not  simply  pushed  later¬ 
ally  by  advancing  projectile  which  can  lead  to  the  defeat  of  the  ar¬ 
mor  by  the  projectile  via  the  so-called  “Wedge  Through  Effect ”  [35]. 

2.4.  Continuum  material  model  implementation  in  a  user-material 
subroutine 

The  unit-cell  based  material  model  described  in  the  previous 
section  is  next  implemented  in  the  material  user  subroutine,  VU- 
MAT,  of  the  commercial  finite-element  program  ABAQUS/Explicit 
[2].  This  subroutine  is  compiled  and  linked  with  the  finite-element 
solver  and  enables  ABAQUS/Explicit  to  obtain  the  needed  informa¬ 
tion  regarding  the  state  of  the  material  and  the  material  mechani¬ 
cal  response  during  each  time  step,  for  each  integration  point  of 
each  element.  In  the  present  work,  first-order  six-node  general- 
purpose  reduced-integration  solid  elements  (ABAQUS/Explicit  des¬ 
ignation  C3D6R)  are  used. 

The  essential  features  of  the  coupling  between  the  ABAQUS/Ex¬ 
plicit  finite-element  solver  and  the  VUMAT  Material  User  Subrou¬ 
tine  at  each  time  increment  at  each  integration  point  of  each 
element  can  be  summarized  as  follows: 

(a)  The  corresponding  previous  time-increment  stresses  and 
material  state  variables  as  well  as  the  current  time-step  deforma¬ 
tion  gradient  are  provided  by  the  ABAQUS/Explicit  finite-element 
solver  to  the  material  subroutine.  In  the  present  work,  the  strain 
components,  two  variables  pertaining  to  the  not-failed/failed  sta¬ 
tus  of  the  filaments  and  one  variable  pertaining  to  the  deletion  sta¬ 
tus  of  the  finite  element  are  used  as  the  state  variables. 

(b)  Using  the  information  provided  in  (a),  and  the  unit-cell 
based  material  model  presented  in  the  previous  section,  the  mate¬ 
rial  stress  state  as  well  as  values  of  the  material  state  variable(s) 
at  the  end  of  the  time  increment  are  determined  within  the  VU¬ 
MAT  and  returned  to  the  ABAQUS/Explicit  finite-element  solver. 
In  addition,  the  changes  in  the  total  internal  and  the  inelastic 
energies  (where  appropriate)  are  computed  and  returned  to  the 
solver. 


Fig.  5.  Typical  final  element  meshes  used  in  the  transient  non-linear  finite-element  analyses  of  the  impact  of  a  full-jacketed  metal  bullet  impact  onto  an  armor-grade 
composite  test  panel.  Due  to  inherit  symmetry  of  the  problem  only  one-quarter  of  the  model  is  analyzed. 
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3.  Validation  of  the  material  model 

As  stated  earlier,  the  main  objective  of  the  present  work  was  to 
develop  a  micro-structure-dependent  physically-based  material 
model  for  the  armor-grade  composites  based  on  0°/90°  cross-plied 
unidirectional  UHMWPE  filaments  and  a  low-content  (<20  mass%) 
polymeric  matrix.  In  this  section,  a  simple  projectile/armor  impact 
problem  is  described.  The  problem  is  used  to  carry  out  preliminary 


Fig.  6.  An  example  of  the  temporal  evaluation  of  the  materials  within  the  22  mm- 
thick  armor  laminate  and  the  M855  projectile.  The  initial  projectile  velocity  equals 
800  m/s. 


testing  and  validation  of  the  proposed  unit-cell  continuum-dam¬ 
age  based  material  model  against  the  previously  reported  experi¬ 
mental  results  of  Iremonger  [36]  who  impacted  five  different 
thickness  (4.2  mm,  11  mm,  15  mm,  22  mm  and  32  mm  with  areal 
densities  given  in  Table  4)  UHMWPE  panels  with  two  types  of  pro¬ 
jectiles.  Iremonger  [36]  fired  two  different  5.56  mm  NATO  projec¬ 
tiles,  the  Ml  93  lead  cored  projectile  with  a  gilding  metal  jacket  and 
the  M855,  which  is  similar  to  the  Ml 93  but  contains  a  hardened 
steel  tip  attached  to  the  lead  core.  The  approximate  masses  of 
the  Ml 93  and  M855  projectiles  are  3.5  g  and  4.0  g,  respectively. 
For  the  present  computational  investigation,  the  panels  were  sim¬ 
ulated  having  a  circular/disk  shape  with  a  radius  of  90  mm. 

An  example  of  the  initial  configuration  of  the  projectile/armor 
finite-element  system  analyzed  here  is  shown  in  Fig.  5.  Due  to 
the  inherent  symmetry  of  the  problem,  only  one-quarter  of  the 
model  is  analyzed  and  the  appropriate  symmetry  boundary  condi¬ 
tions  are  applied  along  the  planes  of  symmetry.  Typically  (one- 
quarter  of)  the  projectile  is  discretized  in  terms  of  1955  first-order 
tetrahedron  elements  C3D4,  while  (one-quarter  of)  the  composite 
panel  is  discretized  using  382  first-order  six-node  brick  elements, 
C3D6,  per  1.1  mm-thick  lamina.  To  reduce  the  computational  bur¬ 
den,  the  size  of  the  C3D6  elements  is  chosen  to  match  that  of  the 
C3D4  elements  only  in  the  region  of  the  composite  panel  impacted 
and  greatly  affected  by  the  bullet.  A  courser  mesh  was  used  in  the 
section  of  the  composite  less  affected  by  the  bullet  impact. 

The  three  metallic  materials  (steel,  lead  and  copper)  present  in 
the  two  types  of  bullets  are  modeled  using  (a  linear  equation  of 
state,  the  Johnson-Cook  strength  model,  the  Johnson-Cook  failure 
model  and  an  erosion  algorithm  based  on  the  maximum  allowable 
instantaneous  geometrical  strain).  Considering  the  fact  that  these 
material  models  were  reviewed  in  our  recent  work  [37],  they  will 
not  be  discussed  any  further  here.  The  panel  was  assigned  the 
material  model  developed  in  the  present  work. 

To  define  the  initial-conditions,  zero  initial  velocities  were  as¬ 
signed  to  all  the  nodes  of  the  composite  panel  while  a  constant 
velocity  in  the  negative  Z  direction  was  assigned  to  all  the  nodes 
of  the  bullet.  Four  initial  projectile  velocities  were  considered: 
600  m/s,  700  m/s,  800  m/s  and  900  m/s.  To  mimic  the  effect  of 
clamping  along  the  armor  edges,  fixed  boundary  conditions  are  ap¬ 
plied  to  all  the  peripheral  nodes  of  the  composite  panel. 

The  same  hard  interaction  algorithm  as  that  used  in  Section  2.1 
was  also  employed  to  model  interactions  between  the  bullet  and 
the  target.  To  account  for  the  sliding-friction  resistance  between 
the  bullet  and  the  target,  a  simple  Coulomb  friction  model  was 
used. 

Computational  analyses  were  run  on  a  machine  with  two 
2.33  GHz  Quad-core  Intel  Xeon  processors  with  16GB  of  RAM.  A 
typical  run  involving  the  11  mm  thick  composite  panel  took 
~12  min  while  in  the  case  of  a  32  mm-thick  panel  the  wall-clock 
computational  time  was  ~45  min. 

4.  Presentation  of  the  results  and  discussion 

As  discussed  earlier,  the  unit-cell  continuum-damage  material 
model  for  0°/90°  cross-ply  UHMWPE  filament  based  armor-grade 
composites  is  validated  by  comparing  the  computational  results 
obtained  in  a  series  of  transient  non-linear  dynamics  finite-ele¬ 
ment  analyses  discussed  in  Section  3  with  their  experimental 
counterparts  obtained  in  the  work  of  Iremonger  [36].  Iremonger 
[36]  investigated  five  armor  panels  with  a  thickness  from  4.2  mm 
to  32  mm,  used  two  types  of  bullets  (M855  and  Ml  93)  and  four  ini¬ 
tial  velocities.  Thus  the  full-factorial  test  matrix  involves 

5.  x  2.  x  4  =  40  experiments.  Not  all  40  experiments  were  con¬ 
ducted  in  the  work  of  Iremonger  [36]  and  only  selected  results 
were  reported  for  the  experiments  that  were  carried  out.  Three 
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types  of  results  were  reported:  (a)  the  success  of  panel  in  stopping 
the  bullet.  Iremonger  [36]  used  the  following  nomenclature:  under¬ 
matched  (to  denote  the  cases  when  the  panel  was  successful  in 
stopping  the  projectile),  overmatched  (to  denote  the  cases  when 
the  panel  was  fully  penetrated  by  the  bullet  but  was  able  to  absorb 
a  substantial  amount  of  the  bullets  kinetic  energy),  and  grossly 
overmatched  (to  denote  the  cases  when  the  bullet  was  able  to  fully 
penetrate  the  panel  by  punching  out  a  circular-disk  shaped  plug  of 
the  panel  material  without  a  significant  loss  in  the  bullets  kinetic 
energy;  (b)  post-mortem  micrographs  of  the  vertical  cut  sections 
of  the  armor  panel  passing  through  the  axis  of  the  penetration 
hole;  and  (c)  temporal-evolution  plots  of  the  back-face  bulge 
height  and  the  bulge  diameter  of  the  impacted  panel.  In  the 
remainder  of  this  section,  a  comparison  is  made  between  the  com¬ 
putational  and  the  experimental  results  for  each  of  these  three 
types  of  results. 

An  example  of  the  typical  computational  results  pertaining  to 
the  temporal  evolution  of  materials  in  the  armor-laminate  and  in 
the  projectiles  is  shown  in  Fig.  6. 

4.1.  The  success  of  armor  test  panels  in  stopping  the  bullet 

The  success  of  different  composite  panels  in  stopping  the  M855 
bullet  at  different  bullet  velocities  as  determined  experimentally 
by  Iremonger  [36]  and  computationally  in  the  present  work  as  dis¬ 


played  in  Table  4.  It  is  seen  that  the  overall  agreement  between  the 
two  sets  of  results  is  quite  reasonable  suggesting  that  the  proposed 
unit-cell  continuum-damage  based  material  model  for  cross-plied 
UHMWPE  filament  based  armor-grade  composites  is  capable  of 
capturing  the  essential  features  of  the  mechanical  behavior  of  this 
material  under  ballistic  loading  conditions.  Results  analogous  to 
those  displayed  in  Table  4  were  not  reported  in  [36]  for  the 
Ml 93  projectile;  hence  a  similar  experiment/computation  compar¬ 
ison  cannot  be  made  for  this  type  of  projectile. 

4.2.  The  mode,  the  extent  and  the  spatial  distribution  of  damage 

In  Ref.  [36],  three  micrographs  of  the  vertical  cut  sections  of  the 
composite  panel  (each  passing  through  the  axis  of  the  penetration 
hole)  were  provided.  Due  to  copyright  restrictions,  only  schematics 
of  these  micrographs  are  included  in  the  present  work.  The  three 
micrographs  correspond  to  the  following  panel/bullet/test  condi¬ 
tions.  (a)  11mm  thick  test  panel/M855/600  m/s,  Fig.  7(a);  (b) 
22  mm  thick  test  panel/M855/800  m/s,  Fig.  8(a);  (c)  22  mm  thick 
test  panel/Ml 93/800  m/s,  Fig.  9(a).  The  corresponding  computa¬ 
tional  counterparts  revealing  the  spatial  distribution  of  damage 
in  the  armor  panel  obtained  in  the  present  work  are  displayed  in 
Figs.  7-9(b),  respectively. 

A  comparison  of  the  results  displayed  in  Figs.  7-9(a)  with  the 
ones  displayed  in  Figs.  7-9(b),  respectively,  reveals  that  the  overall 


Fig.  7.  A  comparison  between  (a  schematic  of)  the  experimental  results  obtained  in  Ref.  [36]  (a),  and  their  computational  counterparts  obtained  in  the  present  work  (b), 
pertaining  to  the  spatial  distribution  of  damage  in  case  of  a  M855  bullet  initially  propelled  at  a  velocity  of  600  m/s  impacting  an  1 1  mm  thick  test  panel. 
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Fig.  8.  A  comparison  between  (a  schematic  of)  the  experimental  results  obtained  in  Ref.  [36]  (a),  and  their  computational  counterparts  obtained  in  the  present  work  (b), 
pertaining  to  the  spatial  distribution  of  damage  in  case  of  a  M855  bullet  initially  propelled  at  a  velocity  of  800  m/s  impacting  a  22  mm  thick  test  panel. 


agreement  between  the  two  sets  of  results  is  reasonable.  More 
specifically: 

(a)  In  the  case  of  1 1  mm-thick  composite  panel  impacted  by  an 
M855  projectile  at  an  initial  velocity  of  600  m/s,  both  the  experi¬ 
ment,  Fig.  7(a),  and  the  numerical  results,  Fig.  7(b),  show  a  “Punch 
through ”  mode  of  penetration  which  is  dominated  by  transverse 
shearing/cutting  of  the  filaments  and  associated  with  relatively 
low  absorption  of  the  projectiles  kinetic  energy.  In  addition,  the 
size  of  penetration  hole  and  its  changes  through  the  armor-panel 
thickness  are  reasonably  well  reproduced  by  the  present  material 
model,  Fig.  7(a)  and  (b). 

(b)  In  the  case  of  a  22  mm-thick  composite  panel  impacted  by 
an  M855  projectile  at  an  initial  velocity  of  800  m/s,  both  the  exper¬ 
iment,  and  the  numerical  results  show  only  a  partial  penetration 
(Fig.  8(a)  and  (b)).  In  the  same  figures  it  is  seen  that  the  depth  of 
the  penetration  hole  is  some  what  under  predicted  by  the  compu¬ 
tational  analyses.  On  the  other  hand,  both  the  experiment  and  the 
computational  analysis  predict  that  initial  penetration  of  the  com¬ 
posite  panel  is  first  dominated  by  filament  shearing/cutting  and 
subsequently  by  filament/matrix  de-bonding/delamination.  It 
should  be  noted  that  the  present  material  model  is  based  on  the 
use  of  a  homogenization  technique  which  yields  an  equivalent  sin¬ 
gle-phase  material  while  the  actual  material  contains  two  phases 
(filaments  and  the  matrix).  Hence,  delamination  in  the  present 
analysis  appears  as  a  region  of  removed  elements  rather  than  a  re¬ 


gion  where  de-bonding  between  adjoining  phases  has  taken  place. 
Despite  these  differences,  the  extent  of  delamination  predicted  by 
the  current  model  appears  comparable  to  that  observed  experi¬ 
mentally.  Also,  the  extents  of  back-face  bulging  observed  experi¬ 
mentally,  Fig.  8(a),  and  the  one  predicted  computationally, 
Fig.  8(b)  are  in  reasonably  good  agreement.  It  should  be  noted  that 
in  the  case  of  the  computational  analysis  some  of  the  elements  at 
the  composite-panel  back-face  in  which  filament  failure  has  taken 
place  attained  a  very  low  level  of  stiffness  which  made  them  un¬ 
dergo  relatively  large  strains.  This  is  the  main  cause  of  the  ob¬ 
served  “bumps”  at  the  armor-panel  back-face. 

(c)  In  the  case  of  a  22  mm-thick  composite  panel  impacted  by 
the  Ml 93  projectile  at  an  initial  velocity  of  800  m/s,  both  the 
experiment,  and  the  numerical  results  show  only  a  partial  penetra¬ 
tion  of  the  target  with  comparable  depth  of  the  penetration  holes, 
Fig.  9(a)  and  (b).  As  in  the  case  of  Fig.  8(a)  and  (b),  the  present  com¬ 
putational  analysis  predicts  reasonably  well  the  extent  of  delami¬ 
nation  within  the  armor  and  the  extent  of  back-face  bulging. 
Again,  low  stiffness  elements  containing  failed  filaments  are  the 
main  cause  of  the  observed  bumps  at  the  armor-panel  back-face, 
Fig.  9(b). 

In  his  work,  Iremonger  [36]  identified  three  distinct  regions  of 
composite  failure  when  subjected  to  impact  from  the  small  caliber 
projectiles:  (a)  an  initial  penetration  phase  dominated  by  fiber 
shearing/cutting  by  the  projectile  and  to  a  larger  extent  by  plastic 
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deformation  of  the  projectile;  (b)  extensive  delamination  of  the 
composite  material  accompanied  by  destabilization  and  break¬ 
up/fragmentation  of  the  projectile;  and  (c)  extensive  stretching/ 
bulging  of  the  armor  back-face  which  was  accompanied  by  exten¬ 
sive  stretching  of  the  fibers  enabling  armor  to  absorb  substantial 
portion  of  the  projectile’s  kinetic  energy.  The  computational  results 
obtained  in  the  present  work  (e.g.  Fig.  8(b))  clearly  confirmed  the 
existence  of  these  three  stages  of  projectile/composite  interaction 
suggesting  that  the  proposed  material  model  for  0°/90°  cross-plied 
UHMWPE  based  armor-grade  composite  is  physically  sound. 

4.3.  Back-face  bulge  topology  and  temporal  evolution 

In  Ref.  [36],  the  temporal  evolution  of  the  composite  panel 
back-face  bulge  height  and  diameter  was  reported  for  only  one 
projectile/bullet/testing  condition.  This  condition  corresponds  to 
a  32  mm  thick  composite  test  panel,  the  M855  bullet  having  an  ini¬ 
tial  velocity  of  900  m/s.  The  results  obtained  in  Ref.  [36]  for  the 
bulge  height  and  for  the  bulge  diameter  are  displayed  in 
Fig.  10(a)  and  (b),  respectively,  and  they  are  denoted  using  filled 
circular  symbols.  The  corresponding  computational  results  ob¬ 
tained  in  the  present  work  are  also  displayed  in  these  figures  and 
they  are  denoted  using  unfilled  circular  symbols.  A  simple  exami¬ 
nation  of  the  results  displayed  in  Fig.  10(a)  and  (b)  reveals  that:  (a) 


while  the  initial  rate  of  increase  of  the  back-face  bulge  height  pre¬ 
dicted  by  the  present  analysis  is  comparable  to  that  measured 
experimentally,  Fig.  10(a),  the  computational  results  under  predict 
the  bulge  height  by  ~5-10  mm.  There  are  several  potential  reasons 
for  this  disagreement,  including:  (i)  the  utilization  of  element  ero¬ 
sion  of  the  projectile  and  target  in  the  current  work;  also  (ii)  the 
extents  of  projectile  damage/fragmentation  were  likely  different 
in  the  experiment  and  in  the  computational  analysis  and  (iii)  both 
the  initial  projectile  velocity  and  the  bulge-height  measurements 
utilized  in  the  work  of  Iremonger  [36]  were  associated  with  exper¬ 
imental  errors  ranging  from  ±5%;  and  (iii).  The  computed  tempo¬ 
ral  evolution  of  the  bulge-height  displayed  in  Fig.  10(a)  shows  a 
decrease  in  the  bulge-height  after  approximately  160  ps.  This  de¬ 
crease  is  associated  with  the  elastic  relaxation  of  the  composite- 
panel  back-face  after  the  projectile  was  defeated  and  pushed  back. 
Similar  observation  was  not  made  by  Iremonger  [36].  Instead,  the 
bulge  height  has  continued  to  increase,  Fig.  10(a).  The  reason  for 
this  discrepancy  is  that  in  the  work  of  Iremonger  [36]  the  projectile 
was  typically  left  buried  within  the  partially  penetrated  armor-pa¬ 
nel  preventing  back-face  elastic  relaxation  and  (b)  except  for  the 
very  initial  stage  of  penetration  by  the  projectile,  the  experimen¬ 
tally  measured  and  computationally  predicted  temporal  evolutions 
of  the  back-face  bulge  diameter  are  in  reasonably  good  agreement, 
Fig.  10(b). 


Fig.  9.  A  comparison  between  (a  schematic  of)  the  experimental  results  obtained  in  Ref.  [36]  (a),  and  their  computational  counterparts  obtained  in  the  present  work  (b), 
pertaining  to  the  spatial  distribution  of  damage  in  case  of  a  Ml 93  bullet  initially  propelled  at  a  velocity  of  800  m/s  impacting  a  22  mm  thick  test  panel. 


M.  Grujicic  et  al./  Composites:  Part  B  40  (2009)  468-482 


479 


Projectile/Armor  Contact  Time,  micro-sec 

Fig.  10.  A  comparison  between  the  experimental  results  obtained  in  Ref.  [36]  and 
their  computational  counterparts  obtained  in  the  present  work  pertaining  to  the 
temporal  evolution  of:  (a)  the  height  and  (b)  the  diameter  of  a  delamination- 
induced  bulge  at  the  composite  back-face  for  the  case  of  a  M855  bullet  initially 
propelled  at  a  velocity  of  900  m/s  impacting  a  32  mm  thick  test  panel. 

4.4.  A  general  discussion  regarding  the  outcome  of  material  model 
validation 

As  stated  earlier,  the  main  objective  of  the  present  work  was  to 
develop,  parameterize  and  validate,  (against  the  relevant  experi¬ 
mental  results),  a  simple  physically-based  computationally-effi- 
cient  continuum-level  material  model  for  a  prototypical  0°/90° 
cross-plied  oriented  polyethylene  fiber-based  armor-grade  com¬ 
posite  material.  While  developing  this  model,  unit-cell  level  fi¬ 
nite-element  analyses  of  the  meso-scale  material  mechanical 
response  and  properties  and  atomic-level  simulations  of  the  fila¬ 
ment/matrix  bonding/de-bonding  had  to  be  employed  giving  the 
present  approach  a  multi-length  scale  character.  The  relevant 
experimental  results  were  taken  from  the  work  of  Iremonger 
[36].  The  material-model  validation  was  carried  out  by  construct¬ 
ing  a  transient  non-linear  dynamics  finite-element  model  consis¬ 


tent  with  the  experimental  setup  used  in  the  work  of  Iremonger 
[36].  The  key  experimental  results  obtained  in  Ref.  [36]  were  then 
compared  with  their  computational  counterparts  obtained  in  the 
present  work  to  judge  the  validity  of  the  present  model.  The  ob¬ 
tained  level  of  qualitative  and  quantitative  agreement  between 
the  two  sets  of  results  suggests  that  the  proposed  material  model 
is  capable  of  capturing  the  essential  behavior  of  a  prototypical 
0°/90°  cross-plied  UHMWPE  based  armor-grade  composite  mate¬ 
rial.  While  the  present  work  was  focusing  on  the  initial  develop¬ 
ment,  parameterization  and  validation  of  the  material  model,  in 
our  future  work,  the  model  will  be  used  to  investigate  in  greater 
details  the  competition  and  interplay  between  various  deforma¬ 
tion,  fracture  and  energy-dissipation  phenomena  which  control 
the  composites  ability  to  defeat  projectiles  by  absorbing  their  ki¬ 
netic  energy. 

It  should  be  noted  that  the  main  projectile  kinetic-energy 
absorption  mechanisms  suggested  by  Iremonger  [36]  and  also  sug¬ 
gested  by  the  present  investigation  are  based  on  filament  cutting 
and  filament/matrix  de-bonding  and  delamination.  As  correctly 
pointed  out  by  one  of  the  reviewers  of  the  present  work,  in  the 
composite  material  under  investigation,  both  the  matrix  and  the 
reinforcement  phases  are  of  polymeric  nature  and  the  role  of  vis¬ 
co-elastic  and  inelastic  (e.g.  plasticity,  crazing,  etc.)  deformation 
modes  must  be  considered.  The  contribution  of  the  latter  energy- 
absorption  mechanisms  was  also  considered  in  the  present  work 
and  it  was  found  to  be  relatively  small.  This  finding  simply  implies 
that,  under  high  loading-rate  conditions  analyzed  in  the  present 
work,  the  response  of  the  constituent  materials  is  dominated  by 
their  elastic  behavior.  While  under  lower  and  moderate  loading 
rates,  visco-elasticity  and  inelasticity  may  play  significant  roles  in 
the  mechanical  response  of  the  material  under  investigation;  these 
roles  are  not  significant  under  ballistic  loading  conditions. 

5.  Work  summary  and  conclusions 

Based  on  the  work  presented  and  discussed  in  the  present  man¬ 
uscript,  the  following  main  summary  remarks  and  conclusions  can 
be  drawn: 

1.  A  simple  unit-cell  continuum-damage  based  material 
model  for  prototypical  cross-plied  unidirectional  UHMWPE  fila¬ 
ment  based  armor-grade  composites  has  been  developed  and 
parameterized. 

2.  To  validate  the  model,  a  series  of  transient  non-linear  dynam¬ 
ics  finite-element  computations  pertaining  to  the  impact  of  either 
an  armor-piercing  (AP)  or  a  non-AP  5.56  mm-caliber  full-metal 
jacketed  bullet  into  a  11-32  mm-thick  armor  test  panel  were  car¬ 
ried  out  and  the  results  obtained  were  compared  with  their  exper¬ 
imental  counterparts  reported  in  Ref.  [36]. 

3.  This  comparison  suggested  that,  for  the  most  part,  the  pro¬ 
posed  model  realistically  accounts  for  the  observed  behavior  of 
the  cross-plied  unidirectional  UHMWPE  filament  based  armor- 
grade  composites  under  the  specified  ballistic  loading  conditions. 
The  good  computation/experiment  agreement  pertains  to  the  suc¬ 
cess  of  the  composite  panels  of  different  areal  densities  in  defeat¬ 
ing  the  bullets  at  different  initial  bullet  velocities,  post-mortem 
spatial  distribution  of  damage  within  the  panel  and  the  temporal 
evolution  of  a  bulge  at  the  back-face  of  the  composite. 

4.  The  computational  analysis  was  also  able  to  clearly  delineate 
three  different  stages  of  composite  penetration  by  the  projectile: 
(a)  initial  filament  shearing/cutting  dominated  stage;  (b)  an  inter¬ 
mediate  stage  characterized  by  pronounced  filament/matrix  de¬ 
bonding/decohesion;  and  (c)  the  final  stage  associated  with  the 
extensive  bulging  of  the  armor-panel  back-face  within  which  pro¬ 
nounced  filament  stretching  leads  to  major  absorption  of  the  pro¬ 
jectile  kinetic  energy. 
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Appendix  A.  Atomic-level  calculations  of  filament/matrix 
bonding/de-bonding 

As  discussed  in  Section  2.1  filament/matrix  bonding/de-bond- 
ing  behavior  and  parameterization  of  this  behavior  are  critical  in¬ 
puts  to  the  unit-cell  based  finite-element  calculations  of  the 
meso-scale  material  mechanical  response.  These  inputs  were  ob¬ 
tained  in  the  present  work  using  an  atomic-level  (i.e.  a  molecu¬ 
lar-statics  based)  computational  analysis.  Some  pertinent  details 
regarding  this  analysis  are  presented  in  this  section.  However, 
since  the  procedure  used  closely  follows  the  one  presented  in 
our  recent  work  [34],  more  details  about  the  present  approach 
can  be  obtained  in  Ref.  [34]. 

Al.  Computational  cell 

An  atomic-level  rectangular  computational  cell  is  first  con¬ 
structed  and  filled  with  aligned  molecules  of  the  polyethylene 
(the  filament  material  and  entangled  molecules  of  urethane  (the 
matrix  material).  One  half  of  the  computational  cell  was  filled  with 
one  material  while  the  other  half  with  the  other,  creating  a  well  de¬ 
fined  filament/matrix  interface.  An  infinite  “crystal”  is  next  con¬ 
structed  by  applying  the  periodic  boundary  conditions  to  all 
three  faces  of  the  computational  cell.  The  cell  axes  (a,  b,  c)  are  ini¬ 
tially  aligned  with  the  three  coordinate  axes  (x,y,z).  The  overall 
initial  orientation  of  the  filament/matrix  interface  was  parallel 
with  the  x—y  plane.  The  construction  of  the  computational  unit  cell 
was  done  using  the  Amorphous  Cell  program  from  Accelrys  [38]. 
The  program  allows  adjustment  in  the  cell  size  to  obtain  the  de¬ 
sired  levels  of  materials  density.  The  atomic  configuration  within 
the  cell  is  then  relaxed/optimized  by  minimizing  its  potential 
energy. 


hyper-surface  on  which  the  atomic  nuclei  move.  In  other  words, 
the  knowledge  of  forcefields  enables  determination  of  the  potential 
energy  of  a  system  in  a  given  configuration.  In  general,  the  poten¬ 
tial  energy  of  a  system  of  interacting  particles  can  be  expressed  as 
a  sum  of  the  valence  (or  bond),  EvaienCe,  cross-term,  Ecwss.term,  and 
non-bond,  Enon_bond ,  interaction  energies  as: 

Etotal  E valence  “F  ECWSs-term  “F  Enon.bond  (Al) 


The  valence  energy  generally  includes  a  bond  stretching  term, 
Eb0nd ,  a  two-bond  angle  term,  Eangie,  a  dihedral  bond-torsion  term, 
Etorsion,  an  inversion  (or  an  out-of-plane  interaction)  term,  Eoop, 
and  a  Urey-Bradlay  term  (involves  interactions  between  two  atoms 
bonded  to  a  common  atom),  Eub,  as: 


- valence 


—  Ebond  +  E( 


angle  ' 


i  +  E0( 


(A2) 


A  schematic  explanation  of  the  first  four  types  of  valence  atomic 
interactions  is  given  in  Fig.  Al. 

The  cross-term  interacting  energy,  Ecr0Ss-term ,  accounts  for  the  ef¬ 
fects  such  as  bond  length  and  angle  changes  caused  by  the  sur¬ 
rounding  atoms  and  generally  includes:  stretch-stretch 
interactions  between  two  adjacent  bonds,  Ebond_bond,  stretch-bend 
interactions  between  a  two-bond  angle  and  one  of  its  bonds, 
Ebond-angie,  bend-bend  interactions  between  two  valence  angles 
associated  with  a  common  vertex  atom,  Eangie-angie ,  stretch-torsion 
interactions  between  a  dihedral  angle  and  one  of  its  end  bonds, 
E  end -bond— torsion*  stretch-torsion  interactions  between  a  dihedral  an¬ 
gle  and  its  middle  bond,  Emiddie_bond_torsiont  bend-torsion  interactions 
between  a  dihedral  angle  and  one  of  its  valence  angles,  Eangie_torsion, 
and  bend-bend-torsion  interactions  between  a  dihedral  angle  and 
its  two  valence  angles,  Eangie_angie-torsion,  terms  as: 


Ebond—bond  “F  Eangie+angieEbond_angie  +  Eend  bond_ torsion 
~F~  Emiddie _bond— torsion  ~F"  Eangie_torsion  +  Eangie_angie_torsion 


(A3) 


The  non-bond  interaction  term,  Enon_bond,  accounts  for  the  inter¬ 
actions  between  non-bonded  atoms  and  includes  the  van  der 
Waals  energy,  Evdw,  the  Coulomb  electrostatic  energy,  EC0uiomb, 
and  the  hydrogen  bond  energy,  EH_bond ,  as: 

Enon-bond  =  Evdw  +  ^Coulomb  “F  ^H-faond  (A4) 


A.2.  Forcefields 

It  is  well-established  that  a  crucial  point  in  the  atomistic  simu¬ 
lations  of  multi-particle  systems  is  the  choice  of  the  forcefields 
which  describe,  in  an  approximate  manner,  the  potential  energy 


Inter-  and  intra-molecular  atomic  interactions  in  the  filament 
and  the  matrix  and  between  the  two  materials  are  modeled,  in 
the  present  work,  using  COMPASS  (Condensed-phased  Optimized 
Molecular  Potential  for  Atomistic  Simulation  Studies),  the  first  ab 
initio  forcefield  that  enables  an  accurate  and  simultaneous  predic- 
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Fig.  Al.  A  schematic  of  the:  (a)  stretch;  (b)  angle;  (c)  torsion;  and  (d)  inversion  valence  atomic  interactions. 
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tion  of  various  gas-phase  and  condensed-phase  properties  of  or¬ 
ganic  and  inorganic  materials  [39,40].  The  specific  functional  forms 
for  Eqs.  (A1)-(A4)  as  defined  within  the  COMPASS  forcefield  can  be 
found  in  Ref.  [34]. 

A3.  Computational  method 

Determination  of  the  filament/matrix  bonding/de-bonding 
behavior  at  the  atomic-level  was  accomplished  in  the  present  work 
by  subjecting  the  computational  cell  to  an  appropriate  deformation 
mode  followed  by  potential-energy  minimization  of  the  cell  with 
respect  to  the  atomic  positions  and  the  free  computational-cell 
edge  lengths  and  angles  (molecular  statics  simulations).  During 
these  simulations,  the  potential  energy  of  the  cell  is  recorded  and 
the  atomic  microstructure  closely  monitored  for  the  onset  of  de¬ 
bonding.  Two  modes  of  deformations  were  used:  (a)  uni-axial  ten¬ 
sion  in  the  direction  normal  to  the  filament/matrix  interface  and 
(b)  simple  shear  in  a  direction  parallel  with  the  filament/matrix 
interface.  All  the  molecular  statics  simulations  are  carried  out 
using  the  Discover  program  by  Accelrys  [41  ]. 

The  potential  energy  minimization  within  Discover  is  carried 
out  by  combining  the  Steepest  Descent,  Conjugate  Gradient  and 
the  Newton’s  minimization  algorithms.  These  algorithms  are  auto¬ 
matically  inactivated/activated  as  the  atomic  configuration  is 
approaching  its  energy  minimum  (i.e.  the  Steepest  Descent  method 
is  activated  at  the  beginning  of  the  energy  minimization  while  the 


Newton’s  method  is  utilized  in  the  later  stages  of  the  simulation). 
Deformation  of  the  computational  cell  was  carried  out  using  a  Dis¬ 
cover  input  file  which  is  written  in  a  BTCL  (Basic  Tool  Command 
Language)  language.  This  enabled  the  use  of  a  scripting  engine  that 
provides  very  precise  control  of  simulation  jobs,  e.g.  a  cell  deforma¬ 
tion  to  be  carried  out  in  small  steps  each  followed  by  energy  min¬ 
imization.  To  help  prevent  the  computational  crystal  from  settling 
into  a  nearby  metastable  higher-energy  configuration,  a  10,000- 
step  298  K  NVT  molecular  dynamics  run  was  introduced  between 
the  cell  deformation  and  the  energy  minimization  steps.  An  exam¬ 
ple  of  the  initial  and  deformed  configurations  of  the  computational 
cell  is  shown  in  Fig.  A2(a)  and  (b). 

The  procedure  described  above  yielded  a  plot  of  the  potential 
energy  increase  (relative  to  that  in  the  initial  optimized  computa¬ 
tional  cell)  vs.  the  extension  Ac  of  the  computational  cell  in  the 
z-direction.  The  gradient  of  the  potential  energy  increase  with  re¬ 
spect  to  the  change  in  the  controlled  computational-cell  parameter 
(i.e.  the  edge  length  or  the  angle)  normalized  by  the  appropriate 
loaded  cross  sectional  area  is  then  used  to  determine  the  corre¬ 
sponding  normal/shear  stress  associated  with  the  given  deforma¬ 
tion  mode.  Likewise,  properly  normalized  changes  in  the  cell 
parameters  are  used  to  compute  the  corresponding  cell  strain(s) 
(i.e.  interfacial  displacement-discontinuities). 

Using  the  aforementioned  procedure,  it  was  possible  to  con¬ 
struct  the  previously-defined  (linear)  filament/matrix  interfacial 
traction  vs.  interfacial  displacement-discontinuity/strain  relations 


Fig.  A2.  An  example  of  the:  (a)  initial  and  (b)  deformed  atomistic  computational  cell  used  in  molecular  statics  analysis  of  the  filament/matrix  bonding/de-bonding  behavior. 
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(one  for  the  normal  and  one  for  the  shear  deformation  modes) 
both  during  elastic  re-loading/unlading  and  loading  associated 
with  interfacial  damage.  For  each  of  the  two  deformation  modes, 
these  relations  are  characterized  by  three  parameters:  (a)  a  criti¬ 
cal  normal  or  shear  interfacial  displacement-discontinuity/strain 
at  which  damage  initiation  begins;  (b)  the  corresponding  normal 
or  shear  traction  (i.e.  the  interfacial  strengths);  and  (c)  a  normal 
or  shear  interfacial  displacement-discontinuity/strain  at  which 
complete  filament/matrix  decohesion  takes  place. 
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